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IMPROVED METHOD AND APPARATUS FOR THE MEASUREMENT OF 

CARDjflC OUTPUT 

The present invention relates to an improved 
5 method and apparatus for the measurement of cardiac 
output and in particular to an improved method and 
apparatus which has a rapid rate of response. 

It has been suggested since at least 1904 
(Erlanger and Hooker, Bull. John Hopkins Hosp. 15:179) 

10 that the arterial pulse pressure could be regarded as 
a rough index to the stroke volume of the heart and, 
in combination with the heart rate, could provide the 
cardiac output. This approach was found to be 
simplistic, particularly at high arterial pressures 

15 where the relationship breaks down. 

Accordingly, various correction factors were 
introduced such as age, weight, height and aortic 
compliance/distensibility . Whilst these correction 
factors improved the results there was still a 

20 disturbing number of conditions in patients where 
disappointing results were obtained, e.g. patients 
with low cardiac outputs and high arterial pulse 
pressures, or vice versa. 

Various approaches to the derivation of cardiac 

25 output from pressure waveform analysis have also been 
attempted: e.g. from the pulse wave velocity and by 
integration of the pulse contour. As stand alone 
methods these have proved equally disappointing with a 
major difficulty in measuring the major morphological 

3 0 features of the waveform in pathologic situations, 
e.g. peaks, systolic areas, dicrotic notch etc. 
During major surgery, or in an intensive care 
situation, the pressure waveform may not exhibit a 
dicrotic notch at all. 

35 In addition to using correction factors Hamilton 
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(The physiology of cardiac output. Circulation 8: 
527, 1953) suggested that cardiac output could be 
derived from a patient's blood pressure pulse height 
following calibration by another cardiac output 
5 method, such as dye dilution. Furthermore, more 
recently Jansen et al (Continuous cardiac output 
monitoring by pulse contour during cardiac surgery, 
Eur. Heart Journal 1990, 11:26-32) discuss a 
•corrected pulse contour 1 method where the pulse 
10 contour is calibrated used a thermodilution technique. 

US Patent No. 540079 3 describes a method for 
determining blood stroke volume from a pulsatile aorta 
blood-flow pressure signal in a subject. The method 
uses a simulation model of the aorta as a transmission 
15 line supplemented with a Windkessel compliance 

including a pressure volume relationship for the aorta 
that is known in the art. 

In essence, the pressure waves recorded from the 
aorta are used to calculate the characteristic 
20 impedance of the transmission line, and the parameters 
of the Windkessel are adapted until the pressure in 
the model is consistent with that recorded in the 
aorta. The flow indicated by the model is then 
integrated over the period of systole. Ideally, this 
25 method requires a high fidelity pressure recording 

system in the aorta. Although a method of correcting 
a pressure measurement in a peripheral artery is 
mentioned, this method cannot be used with the poor 
frequency response given by most pressure transducers 
30 that are now routinely in clinical use: in the 

presence of noise an "anti-resonance f ilter" cannot 
recover the information that is lost by the poor 
quality of these transducers. 

Finally, Irlbeck et al (Continuous measurement of 
35 cardiac output with pulse contour analysis, 
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Anaesthetist 1995, 44:493-500) propose that the 
calculation of stroke volume and cardiac output is 
possible by analysing the area under the systolic part 
of the arterial pulse pressure waveform, together with 
5 an individual calibration factor (Zao) to account for 
individually variable vascular impedance. 

Unfortunately, this analysis is of limited 
applicability to clinical measurement since the 
waveforms recorded for the patients for whom it would 

10 be really useful often do not show the features 

reguired for the recognition of the systolic part of 
the pressure waveform. The pulse is often irregular 
and each heartbeat may generate several minor waves 
superimposed upon the primary pulse waveform. 

15 W e have now developed an improved method and 

apparatus for measuring cardiac output which does not 
suffer the aforementioned disadvantages of the prior 
art methods discussed above. 

Accordingly, in a first aspect the present 

20 invention provides a method for the measurement of 
cardiac output in a patient, which method comprises 
the steps of:- 

(i) recording and storing the arterial blood 
pressure waveform of a patient from a 

25 blood pressure monitoring device over a 

period of time; 

(ii) subjecting the waveform obtained in step 
(i) to a non-linear transformation that 
corrects for the variation of the 

30 characteristics of the arterial system 

with pressure; 

(iii) subjecting the corrected waveform curve 
from step (ii) to autocorrelation in order 
to derive the pulsatility and heart rate 

35 of the corrected waveform; 
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(iv) calculating the nominal stroke volume from 
the pulsatility; and 

(v) obtaining the nominal cardiac output by 
multiplying the nominal stroke volume by 

5 the heart rate. 

In a second aspect the present invention provides 
a method for the measurement of cardiac output in a 
patient, which method comprises the steps of:- 

(a) recording and storing the arterial blood 
pressure waveform of a patient from a 
blood pressure monitoring device over a 
period of time; 

(b) subtracting the mean of the waveform from 
step (a) and subjecting the data so 

15 obtained to autocorrelation; 

(c) transforming the data from step (b) into 
data which relates to the pulsatility and 
heart rate of the waveform; 

(d) calculating the nominal stroke volume from 
20 the pulsatility; and 

(e) obtaining the nominal cardiac output by 
multiplying the nominal stroke volume by 
the heart rate. 

In carrying out either of the methods of the 
25 present invention the arterial blood pressure of a 
patient is monitored continuously by conventional 
means via an arterial line and pressure transducer, 
for example using a Hewlett Packard blood pressure 
monitoring machine. The output data from the machine 
30 is amplified and smoothed to provide a continuous 
blood pressure waveform. The continuous blood 
pressure waveform is monitored over a period of time, 
for example at least two patient heartbeats and 
generally for a period of ten seconds or less, 
35 preferably for a period of up to four seconds. 
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In carrying out the first method of the present 
invention the pressure is subjected to a non-linear 
transformation before the autocorrelation step whereas 
in other methods correction factors are generally 
5 applied at the end of the calculations. 

As mentioned above, the pressure waveform 
obtained in step (i) of the first method is 
transformed, preferably via a 'look up' curve, with the 
mean of the' data then being found and subtracted, into 

10 data which represents the pressure /volume relationship 
of the arterial system. The basic approximation to a 
look up table is known in the art and the relationship 
is non-linear, a series of such curves being described 
in Remington et al, 1948, Am. J. Physiol 153: 298-308: 

15 Volume elasticity characteristics of the human aorta 
and prediction of the stroke volume from the pressure 
pulse. 

The corrected waveform from step (ii) is then 
subjected to autocorrelation in order to determine the 

20 pulsatility and heart rate of the transformed 

waveform. Autocorrelation is defined in Dictionary of 
Science and Technology, Academic Press, 1992, p 186. 

Autocorrelation is performed by a process of 
multiplication and summation. Firstly, the mean of 

25 the series is found and subtracted from each of the 

values so that the series represents the signal moving 
positive and negative around a mean of zero. 

Tau (t) is used to denote a time shift so that at 
Tau = 0 each of the values of y is multiplied by 

30 itself and all of the products added together. When 
Tau = 1 the first value is multiplied by the second, 
the second by the third, and so on down the series 
representing the record. The values obtained for each 
Tau are then divided by the number of points used in 

35 the correlation (say n/2) and the correlation function 
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plotted as a function of Tau. 

When Tau = 0 the function has a maximum value 
since each of the values is being squared and the 
negative values become positive. If the original 
5 series represented white noise, any non-concurrent 
points would have a random relationship and be as 
likely to be positive as negative. When Tau is 
greater than zero the sum of the products would 
therefore be close to zerof However, if there is any 

10 systematic waveform contained in the original series 
then at an appropriate value for Tau, like is again 
being multiplied by like, and the autocorrelation 
rises to a peak. 

This form of correlation is very powerful in 

15 detecting periodicities; the Tau value of the first 
correlation peak represents the wavelength of the 
periodicity. Depending upon the number of cycles 
sampled by the original series, further peaks may be 
seen representing the Tau shift at which the first 

20 waveform is being multiplied by the next but one, and 
so on- Therefore the interval for one heart beat can 
be found using autocorrelation - it is equal to the 
smallest value of Tau that gives a peak in the 
autocorrelation which is comparable to the mean square 

25 value at Tau - O. 

The value of the function of Tau = 0 is the mean 
square of the original record. Thus if the original 
trace were a series of cycles of a sine wave, the 
autocorrelation would resemble a cosine waveform, the 

30 peak values of which would be the mean square of the 
original sine wave. Taking the square root of the 
peak value gives the root mean square of the original 
sine wave. 

The autocorrelation (of the transformed data with 
35 the mean subtracted) is also used to find the 



WO 97/24982 



PCT/GB97/00017 



- 7 - 



pulsatility. We define pulsatility as 

pulsatility = ylRiS) * f%x~T) 

5 where R(t) is the autocorrelation function and r nin is 
the value of t at which the autocorrelation is at a 
minimum. 

Nominal stroke volume is then found by 
multiplying the pulsatility by a constant. The 

10 nominal cardiac output is obtained by multiplying the 
nominal stroke volume by the heart rate. It will be 
understood that the nominal stroke volume and the 
nominal cardiac output are uncalibrated and may be 
converted into calibrated data, if desired. 

15 The stroke volume is calculated in step (iv) of 

the method , generally by multiplying the pulsatility 
by a calibration factor obtained for a particular 
patient using a known accurate calibration method, for 
example using an indicator dilution method or 

20 thermodilution method. An indicator dilution method 
is described, for example, in W093/09427. The method 
as described in W093/09427 is highly repeatable and 
only one single point calibration is required to give 
the cardiac output. It will be understood, however, 

25 that the method of the present invention may be used 
without calibration in order to show trends in or 
directions of change of the cardiac output of a 
patient. 

Once the adjustment/calibration has been made 
3 0 then the cardiac output can be obtained by multiplying 
the stroke volume by the heart rate. The cardiac 
output can be provided as an appropriate value of few 
seconds coverage and as data, thereby giving a close 
to realtime analysis of rapidly changing events. For 
35 example, the method of the present invention can be 



WO 97/24982 



- 8 - 



PCT/GB97/00017 



used to monitor the cardiac output change following 
the administration of fluids, or to set a temporary 
pacemaker to an optimal rate, or to determine when 
administration of a vasoactive drug may be required. 
5 Another method of adjusting /calibrating the 

pulsatility is by the thermodilution method but this 
method is not preferred since it requires multiple 
readings viz the averaging of at least three points in 
order to give a single stroke volume from which the 
10 cardiac output can be calculated and presents a 
significant risk to the patient. 

In carrying out the second method of the present 
invention steps (a) , (d) and (e) correspond generally 
to steps (i), (iv) and (v) of the first method as 
15 discussed above. Following step (a), the mean of the 
waveform is subtracted and the data so obtained is 
then subjected to autocorrelation in the manner as 
discussed above. The autocorrelated data then 
transforms into data which relates to the pulsatility 
and heart rate of the waveform. The calculation of 
the pulsatility using this method will require scaling 
according to the mean blood pressure, the scaling 
factor being equal to the slope of the pressure-volume 
relationship at the mean pressure. 
25 Although the methods of the present invention 

have been described above generally assuring that a 
plurality of heart beats will be measured, it is 
possible to process single heart beats if they can be 
identified separately. Methods for identifying single 
heart beats are known in the art. In order to find the 
stroke volume of a single beat the method would be 
essentially the same as that described above, but to 
carry out the autocorrelation it would be assumed that 
the blood pressure repeated itself after the end of 
35 the beat. 



20 



30 
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In a further aspect the present invention 
provides a first apparatus for the measurement of 
cardiac output in a patient, which comprises: 

means for monitoring the arterial blood pressure 
5 of a patient; 

means for monitoring the heart rate of a patient; 

means for recording and storing the arterial 

blood pressure waveform of a patient; 

means for transforming the pressure waveform to 
10 correct for the variation of the arterial system 

with pressure; 

means for autocorrelating the transformed data 

and deriving the pulsatility and heart rate of 

the waveform therefrom; and 
15 means to calculate the nominal cardiac output 

from the measurements of heart rate of the 

patient and the pulsatility of the waveform. 

In the apparatus of the invention the means for 
transforming the pressure waveform preferably 
20 comprises a look up table as described above. The 
apparatus may also include means to calibrate the 
device, for example using the indicator dilution 
methods as discussed above. 

In a still further aspect the present invention 
25 provides a second apparatus for the measurement of 
cardiac output in a patient, which comprises: 

means for monitoring the heart rate of a patient; 

means for recording and storing the arterial 

blood pressure waveform of a patient; 
30 means for subtracting the mean of the blood 

pressure waveform and subjecting the data to 

autocorrelation; 

means for transforming the autocorrelated data 
into data relating to the pulsatility and heart 
35 rate of the waveform; and 
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means to calculate the nominal cardiac output 
from the measurements of heart rate of the 
patient and the pulsatility of the waveform. 
Apparatus for carrying out the present invention 
5 may comprise any suitably programmed computer such as 
an IBM compatible computer or a Macintosh personal 
computer. For the analysis of indicator dilution 
curves, the apparatus for measuring cardiac output 
will include a blood sensor such as that disclosed in 
10 WO93/09427. The output of the sensor must then be 
digitized by an analogue-to-digital converter and 
input into the computer. For any other form of 
signals representing physical parameters, similarly, 
the output of some sort of transducer must be 
15 digitized to provide signals for analysis by the 
software loaded in a computer. 

The computer program running on the computer may 
either display the results or can output this to some 
other device. 

20 The present invention will be further described 

with reference to the accompanying drawings in which: 

Figures 1 to 4 show the results obtained using 
method and apparatus of the present invention to 
analyse the arterial blood pressure recordings from 

25 four patients; and 

Figures 5a and 5b compare the cardiac output of 
two patients measured using the method and apparatus 
of the present invention compared to the prior art 
method as disclosed in WO 93/09427. 

30 Referring to the drawings, in each of Figures 1 

to 4 the top left window shows four seconds of 
arterial blood pressure recorded from a patient as the 
upper trace and below this the autocorrelation. In 
the upper right window the heart rate (HR) , cardiac 

35 output (CO) and -pulsatility" (PP) are displayed. The 
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bottom window shows a right to left scrolling record 
of previous history of the cardiac output. 

Fig. - 1 The arterial wave form is smooth and 
regular. The heart rate is 7 9.9 beats per minute. 
5 Fig. 2 - The arterial wave form shows a secondary 

wave at the peak as well as a dicrotic notch on the 
descending limb. The heart rate again is 79.9 beats 
per minute. Both these patients were electrically 
— paced at a nominal 80 bpm. 

10 Fig. 3 - This shows a pronounced and delayed 

diastolic wave. The autocorrelation gives the correct 
heart rate. 

Fig. 4 - This shows an alternating pattern in the 
heart beat. The autocorrelation estimates the rate 

15 correctly and shows that alternate beats correlate 
better than adjacent beats. 

Figures 5a and 5b are graphs which compare 
cardiac output measurements obtained using the prior 
art lithium dilution method as disclosed in WO 

20 93/09427 and the first method of the present 

invention. The close correlation of the results 
clearly indicates the merits of the present invention. 



WO 97/24982 



- 12 - 



PCT/GB97/00O17 



CLAIMS: 

1. A method for the measurement of cardiac 
output in a patient, which method comprises the steps 
5 of:- 

(i) recording and storing the arterial 

blood pressure waveform of a patient 
from a blood pressure monitoring device 

"~~ over a period of time; 

10 (iij subjecting the waveform obtained in 

step (i) to a non-linear transformation 
that corrects for the variation of the 
characteristics of the arterial system 
with pressure; 

15 subjecting the corrected waveform from 

step (ii) to autocorrelation in order 
to derive the pulsatility and heart 
rate of the corrected waveform; 

(iv) calculating the nominal stroke volume 
20 from the pulsatility; and 

(v) obtaining the nominal cardiac output by 
multiplying the nominal stroke volume 
by the heart rate. 

25 2. A method as claimed in claim 1 wherein the 

arterial blood pressure is plotted in step (i) for a 
period of up to ten seconds. 

3. a method as claimed in claim 2 wherein the 
30 arterial blood pressure is plotted in step (i) for a 

period of up to four seconds. 

4. A method as claimed in any one of the 
preceding claims wherein the transformation in step 

35 (ii) is effected using a look up table, with the mean 
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of the data then being found and subtracted. 

5. A method as claimed in any one of the 
preceding claims wherein the nominal stroke volume is 
obtained by multiplying the pulsatility by a 
calibration factor. 

6. A method as claimed in claim 5 wherein the 
calibration obtained for a particular patient using a 
known calibration method. 

7. A method as claimed in claim 6 wherein the 
calibration method is an indicator dilution method. 

5 8. A method for the measurement of cardiac 

output in a patient, which method comprises the steps 
of:- 

(a) recording and storing the arterial 
blood pressure waveform of a patient 

0 from a blood pressure monitoring device 

over a period of time; 

(b) subtracting the mean of the waveform 
from step (a) and subjecting the data 
so obtained to autocorrelation; 

5 ( C ) transforming the data from step (b) 

into data which relates to the 
pulsatility and heart rate of the 
waveform; 

(d) calculating the nominal stroke volume 
0 from the pulsatility; and 

(e) obtaining the nominal cardiac output by 
multiplying the nominal stroke volume 
by the heart rate. 



5 



9. A method as claimed in claim 8 wherein the 
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arterial blood pressure is plotted in step (a) for a 
period of up to ten seconds. 

10. A method as claimed in claim 9 wherein the 
5 arterial blood pressure is plotted in step (a) for a 

period of up to four seconds. 

11. A method as claimed in any one of claims 8 to 
10 wherein the nominal stroke volume is obtained by 

10 multiplying the pulsatility by a calibration factor. 

12. A method as claimed in claim 11 wherein the 
calibration factor is obtained for a particular 
patient using a known calibration method. 



15 



13. A method as claimed in claim 12 wherein the 
calibration method is an indicator dilution method. 

14. Apparatus for the measurement of cardiac 
20 output in a patient, which comprises: 

means for monitoring the arterial blood pressure 
of a patient; 

means for monitoring the heart rate of a patient; 

means for recording and storing the arterial 
2 5 blood pressure waveform of a patient; 

means for transforming the pressure waveform to 

correct for the variation of the characteristics 

of the arterial system with pressure; 

means for autocorrelating the transformed data 
30 and deriving the pulsatility and heart rate of 

the waveform therefrom; and 

means to calculate the nominal cardiac output 
from the measurements of heart rate of the 
patient and the pulsatility of the waveform. 



35 
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15. Apparatus as claimed in claim 14 wherein the 
means for transforming the pressure waveform comprises 
a look up table. 

5 16. Apparatus as claimed in claim 14 or claim 15 

which includes means to calibrate the device. 

17 . Apparatus for the measurement of cardiac 
output in a patient, which comprises: 
10 means for monitoring the arterial blood pressure 

of a patient; 

means for monitoring the heart rate of a patient; 
means for recording and storing the arterial 
blood pressure waveform of a patient; 
!5 means for subtracting the mean of the blood 

pressure waveform and subjecting the data to 
autocorrelation ; 

means for transforming the autocorrelated data 
into data relating to the pulsatility and heart 
20 rate of the waveform; and 

means to calculate the nominal cardiac output 
from the measurements of heart rate of the 
patient and the pulsatility of the waveform. 
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